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ABSTElACr 


Analytical  ej<pressions  for  the  horizontal  and  vertical  forces 
exerted  by  a  floating  ice  sheet  against  a  tvwD-dimensional  sloping  struc¬ 
ture  are  derived  under  the  assunption  that  ice  is  an  elastic,  homogeneous 
material.  The  specific  conditions  investigated  are 

1.  Quasi-steady  form  of  the  governing  equation. 

a.  The  free  end  of  the  seni-infinite  ice  sheet  is 
subjected  to  a  horizontal  force,  a  vertical  force, 
and  a  bending  rroient.  Failure  of  the  ice  sheet 
occurs  before  its  free  edge  becomes  fully  energed 
from  or  submerged  into  the  water. 

b.  The  free  end  of  a  s^ni-infinite  ice  sheet  is  sub¬ 
jected  only  to  a  vertical  force,  but  failure  occurs 
after  the  free  end  has  fully  emerged  from  or  submerged 
into  the  water. 

c.  Finite  ice  sheet  subjected  only  to  a  vertical  force 
at  one  free  end  wiiich  remains  at  the  water  surface 
until  failtore  occurs. 

2.  Unsteady  form  of  the  governing  equation  for  a  seid-infinite 
sheet  subjected  to  only  a  vertical  force  at  its  free  end  vdiich  remains  at 
the  water  surface  until  failure  occurs. 

It  appears  that,  for  practical  applications,  the  effect  of  the 
horizontal  force  and  bending  mcanent  at  the  free  end  of  the  ice  sheet  on 
the  sheet  deformation  is  negligible.  It  is  shown  that  a  finite  ice  sheet 
of  length  at  least  three  times  a  characteristic  length  e?q)ressed  in  terms 
of  the  ice  mechanical  properties  can  be  treated  as  a  sani-infinite  ice 
sheet.  Also  it  is  shown  that  when  the  ice  thickness  is  smaller  than  a 
critical  value  which  depends  upon  the  ice  mechanical  properties  and  the 
direction  (upward  or  downward)  of  the  vertical  force  exerted  at  the  sheet 
free  end,  the  failure  force  of  the  ice  sheet  is  independent  of  the  ice 
elastic  medulus. 
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The  analytical  results  have  been  verified  experimentally,  and 
used  in  the  experimental  deteirnination  of  the  bending  strength  and  strain 
modulus  of  urea“doped  ice. 
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I .  INTRODUCTION 


HJasi  wade  structures  erected  in  the  high  latitude  regions  of 
the  globe  must  be  built  so  as  to  withstand  the  large  forces  exerted  by 
the  ice  ■which  they  may  come  in  contact  with.  These  structures  may  be 
stationary,  bridge  piers,  warves,  dykes,  light  houses,  drilling  islands 
etc. ,  or  mobile,  tugs,  icebreakers  etc.  The  ice  these  structures  may 
have  to  interact  with  can  be  monolithic  sheets  or  broken  ice  masses 
either  unconsolidated,  ice  rubble,  or  partially  or  totally  refrozen,  ice 
ridges.  Since  the  mechanical  properties  of  monolithic  ice  show  its 
bending  strength  to  be  about  half  its  crushing  strength  (e.g. ,  Pounder 
1965),  ice  thrust  on  structures  can  be  reduced  by  inclining  the  face  of 
the  structiire  to  force  the  ice  to  fail  by  bending. 

Because  of  its  importance  it  is  not  surprising  that  numerous 
studies  have  been  denoted  to  the  problem  of  ice  forces  on  inclined  struc¬ 
tures,  and  related  topics  by  an  analytical  approach,  experimental  simula¬ 
tion  and  combinations  of  both.  For  exaxrple,  Nevel  (1979)  and  Sorensen 
(1978)  have  followed  the  analytical  route,  wiiile  Tryde  (1972),  Ralston  (1979) 
Frederking  (1979),  among  others,  have  followed  an  enpirio-analytical 
approach.  Iifost  of  these  researchers  have  treated  the  case  of  a  three- 
dimensional  structure  interacting  with  an  infinite  floating  ice  sheet, 

■vdiere  three-dimensional  effects  occur,  and  for  which  the  analysis  becomes 
extremely  complicated,  even  ■when  only  quasi  steady  conditions  are  considered. 
The  problem  treated  herein  is  tliat  of  a  free-floating  ice  sheet  pushing 
against  a  two-dimensional  inclined  plane.  The  problem  is  tractable 
mathonatically,  and  it  is  believed  that  the  results  obtained  can  give 
better  insight  into  the  general  phenomenon  of  ice  forces  on  structiures  of 
more  complicated  shape.  In  particular  it  is  shown  that  the  relationship 
between  the  maximum  force  exerted  by  an  ice  sheet  on  an  inclined  plane  and 
the  ice  properties,  bending  strength  and  elastic  modulus,  depends  on 
\diether  the  ice  sheet  is  lifted  out  of  the  ■water  or  pushed  down  into  it, 
or  ■whether  failure  of  the  ice  sheet  occurs  before  or  after  its  frontal 
edge  beccmes  fully  anerged  from  or  submerged  into  the  water,  and  also 
on  how  high  the  frontal  edge  emerges  or  how  deep  it  submerges  before 
failure  occurs. 
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II. 


STA1 


T  OF  PROBLBvl 


The  problem  to  be  considered  is  that  of  an  ice  sheet  floating 
at  the  surface  of  a  body  of  water  and  which,  under  the  influence  of 
external  forces,  such  as  wind  shear,  is  climbing  up  or  down  a  sloping 
stmcture.  Only  the  two-dimensional  problem  is  considered  here.  The 
ice  sheet  of  uniform  thickness  h,  and  width  b  is  assumed  to  be  hcxnogeneovis 
and  to  behave  as  an  elastic  material  characterized  by  a  modulus  of 
elasticity  E,  a  bending  strength  and  a  ccxipression  strength  o^. 

The  chosen  systan  of  coordinates  has  its  origin  on  the  line  of 
intersection  between  the  water  surface  and  the  end  of  the  ice  sheet. 

The  x-axis  is  horizontal  in  the  plane  of  the  water  surface  in  the  opposite 
direction  of  motion  of  the  ice;  the  y-axis  is  vertical  upward  (see 
Figure  1).  The  origin  of  time  is  chosen  as  the  instant  when  the  ice 
strikes  the  inclined  structure.  The  horizontal  di^lacement  velocity  of 
the  ice,  V  ,  will  be  assumed  to  ranain  constant  throughout  the  process. 

For  time  t  _<  0,  the  buoyancy  force  balances  the  ice  weight, 
such  that  the  iimersed  depth  of  ice,  h^,  is  given  by 

h  =  —  h  (1) 

O  D 
^W 

where  and  are  the  densities  of  ice  and  water,  respectively. 

For  time  t  >  0,  as  the  ice  climbs  up  the  structure,  it  is 
subjected  to  a  reaction  force  R  of  coiponents  N  and  T  in  the  normal  and 
tangential  directions,  respectively,  relative  to  the  structure,  and  of 
ccarponents  P  and  F  in  the  xy-systan  of  coordinates  (see  Figure  1),  such 
that 


F  =  N  cos  e  -  T  sin  0  (2) 

P  =  N  sin  e  +  T  cos  0 

where  0  is  the  angle  of  inclination  of  the  stmcture  mth  respect  to  the 
horizontal  plane. 
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Figure  1.  General  Definition  Sketch 
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When  the  dynamic  friction  coefficient,  y,  at  the  ice-structure 
interface  is  introduced  then 


T  =  pN 

and 

p  =  sy9_+  ycose  Y  =  <i,F 

cose  -  ysine  ^ 


(3) 


The  reaction  force  R  acts  at  the  point  of  contact  of  the  ice  sheet  with 
the  structure. 

Under  the  action  of  F  and  P,  the  ice  sheet  xindergoes  a  dis¬ 
placement  y(x,t)  frcxn  its  initial  free  floating  position,  and  is  sub¬ 
jected  to  a  moment  distribution  M(x,t).  The  ice  sheet  will  fail  at  a 
distance  x^  from  the  upstream  end,  whe3?e  the  stress  at  the  top  or  botton 
of  the  ice  sheet  reaches  either  the  bending  strength,  or  the  ccxn- 

pressive  strength,  a  ,  of  the  ice.  However,  since  the  tensile  strength 

c 

of  ice  is  approximately  half  its  coipressive  strength  (see,  for  exanple. 
Pounder  1965,  Weeks  and  Assur  1969),  it  will  be  assumed  that  the  ice 
sheet  will  fail  in  tension.  The  distance  of  failure  x^  is  then  determined 

by 


_Mfh  p 

°b  21  bh 


(4) 


or 

=  ^  _  iZ 

%  21  bh 


(4a) 


bh" 


where  =  M(x^,t)  and  I  =  is  the  nonent  of  inertia  of  the  beam 
cross-section  about  the  neutral  axis. 
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III.  PROBLEM  ANALYSIS 


A.  Governing  Eqviations 

I.  Ice  Sheet  in  Contact  with  Water  Along  its  Whole  Length 

Consider  an  elenentary  slice  of  the  ice  sheet  of  length  dx  as 
shown  in  figure  2a.  The  ice  element  is  displaced  by  the  distance  y 
from  its  free  floating  position,  and  is  subjected  to  the  various  forces 
and  manents  indicated  on  figure  2a.  Application  of  Newton's  second  law 
in  the  y-direction  leads  to 

2 

(p  .bhdx)  ^  =  f  +  dB  -  dW  -  (f  +  dx)  (5) 

at"' 

vdiere  f  =  shear  force 

dB  =  buoyancy  force  on  ice  elanent  =  y  (h  -y)bdx 

wo 

dW  =  weight  of  ice  elonent  =  y  .  hdbx  =  v  tdx 

'  1  'wo 


Equation  5  gives 


bh  (?-  +  i  H) 


(6) 


The  shear  force  distribution  along  the  ice  sheet  is  then  given  by 

f(x,t)  =  F  +  /||dx  (7) 

The  transverse  deformation  y(x,t)  of  the  ice  sheet  is  related  to  the 
nonent  M(x,t)  exerted  at  any  section  by 


M(x,t)  _  1  _ 
El  r 


^  3x 


(8) 


vdaere  r  is  the  radius  of  curvature  of  the  neutral  axis.  When  it  is 
assumed  that  the  transverse  deformation  is  amll  enough  for  (9y/9x)" 
to  be  negligible  with  respect  to  1,  equation  8  reduces  to 
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M  _  3  y 

El  ,  2 
3x 


(9) 


The  reactions  F  and  P  exerted  at  the  bottom  of  the  free  end  in  contact 
with  the  structirre  can  be  replaced  by  the  same  two  forces  exerted  at  the 
neutral  axis  of  the  ice  sheet  plios  a  bending  nonent  as  shown  in 
Figure  1.  When  mcment  equilibriijm  is  applied  to  a  sheet  elanent  of 
length  dx,  or  to  a  segment  of  ice  sheet  extending  fixsn  the  free  end  to  a 
distance  x,  the  following  relationship  for  the  bending  racxnent  exerted  at 
any  section  x  of  the  ice  sheet  is  obtained 


X 

M(x,t)  =  M  +  P(y  -y)  +  /  f  dx  (10) 

o  o  0 

where  y  is  the  vertical  displacement  of  the  free  end  firon  its  free 
floating  position.  When  equation  10  is  differentiated  twice  with  respect 
to  X,  it  yields 


(11) 


When  equations  4,  5  and  9  are  substituted  into  equation  11,  the  governing 
equation  for  the  deformation  y(x,t)  is  obtained  as 


ifz  +  iZ  +IzL 

3X^  ^^3x2  El 


(y  + 


h  .2 

_0  =  0 
g  2^  ^ 

^  3t 


The  following  dimensionless  variables  are  introduced 
-  =  y/L,  Z  =  x/L,  t  =  t/T,  F  =  F/r^bL^ 

M  =  M/y  bL^ 
w 


(12) 


where  L  is  a  characteristic  length  defined  by 


=  (— ) 
V  b-’ 

'  w 


=  (i?-) 


3y. 


w 


(13) 


and  T  is  a  characteristic  time  defined  by 
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T  =  (h^/g)^  (14) 


Then  the  dimensionless  governing  equation  is; 


4r-  2-  2- 

9y  +  4j,F  +  4(y  +  IJ)  =  0  (15) 

3X  3x  3t 


2,  Partially  Emerged  Ice  Sheet 

The  above  derivation  and  the  resulting  governing  equation  15 
are  valid  only  as  long  as  the  ice  sheet  at  any  point  along  its  length  is 
neither  fully  anerged  fron  nor  fiolly  subraerged  into  the  water,  i.e.,  for 


-(h-h^)  <  y(x,t)  <  h^ 


In  particialar,  if  the  free  end  anerges  fron  the  water  prior  to  failure, 
the  above  analysis  needs  to  be  modified  over  the  emerged  length  x^  as 
shown  in  Figure  2b.  Newton’s  second  law  applied  to  an  elementary  ice 
elanent  dx  in  the  region  x  <  x^  yields 

2 

(p^bh  dx)  =  f  -  dW  -  (f  +  1^)  dx  (5a) 


or 


(6a) 


Equations  10  and  11  ranain  valid  over  0  <  x  <  x^,  therefore  the  governing 
equation  for  the  displacanent  y(x,t)  over  the  anerged  length  of  ice 
sheet  becomes 


=  0 


(12a) 


or,  in  dimensionless  form. 


^  +  44.F  -^  +  4(h  +  i^)  =  0 
3X  3x  3t 


^\bere  h  =  h^/L  and,  valid  for  0  <  x  <  x 


o 


(15a) 
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Eqxiations  12  and  15  ramin  valid  for  x  >  x^,  i.e.,  x  >  x^. 


3.  Final  Ranarks 

In  the  derivation  of  Equations  15  and  15a  it  was  inplicitly 
assumed  that  the  vertical  force  F  exerted  at  the  free  end  of  the  ice 
sheet  was  directed  upward.  If,  on  the  contrary,  F  is  directed  downward, 
these  equations  are  still  valid  provided  that  the  displacanent  y 
(x,t)  be  taken  positive  downward  and  that  h^  be  replaced  by  (^-5^)  in 
Equation  15a.  In  sunrnary,  the  dimensionless  governing  eqioations  are: 


+  4  +  4(y  +  2^)  =  0 

9x  ax  at 


for  y  <  H 


^  +  4  <f)F  ^  +  4(H  +  =  0 

ax  ax  at 


for  y  >  H 


(15) 


(16) 


^^here  H  =  h^  if  ice  sheet  subjected  to  a  lift  force. 


H  =  h-h^  if  ice  sheet  subjected  to  a  downward  push 


B.  Initial  and  Boijndary  Conditions 

1.  Initial  Conditions 

When  the  unsteady  form  of  the  governing  equations  are  considered, 
the  initial  condition  is  sinply  y(x,0)  =  0. 

2.  Boundary  Conditions 

There  are  three  possible  boundaries  to  be  considered;  the 

upstream  edge  of  the  ice  sheet,  x  =  0,  which  is  in  contact  with  the 

structure  and  is  submitted  to  the  forces  F  and  P,  and  a  nonent  M  ; 

o 

the  downstream  end  of  the  sheet  corresponding  to  x  =  °°  for  a  sarii- 
infinite  sheet,  and  to  x  =  for  a  sheet  of  finite  length  and  finally, 
in  the  case  of  a  partially  emerged  (or  submerged)  sheet,  the  location  x 

=  X 

o 
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a)  at  X  =  0  the  boundary  conditions  are 


y  =  t 


li 

3x2 

rN  3— 

9  y 
35" 


4  M 


+  4  (j)F 


■p  iZ  =  4  F 


9x 


(17) 

(18) 

(19) 


\diere  is  the  dimensionless  constant  vertical  velocity  of  motion  of 
ice  edge  in  contact  with  the  structure 


b)  at  X  ~  for  a  sani-infinite  ice  sheet 


3n- 

— ^  ^  0  n  =  0,1,2,... 

9x 

c)  at  X  =  t  for  a  sheet  of  finite  length 

=  0 


2— 

M(x  =  t)  =  0  or 

9x2 1 


x=)l 


(20) 


(21) 


d)  at  X  =  x^,  for  the  case  of  an  ice  sheet  totally  emerged,  or  submerged, 
over  a  length  x^.  Let  y^(x,t)  be  the  solution  of  Equation  (16),  valid 
for  X  <  x^,  and  ^^(x,!)  the  solution  of  Equation  15,  valid  for  x  >  x^. 

We  shall  require  continuity  of  displacement,  slope,  moment  and  shear,  at 
X  =  x^  namely 


y 


1 


y 


=  H 


(22a) 


!Zi  =!!2 

9x  9x 


9x2 


(22b) 


(22c) 
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-3  -3 

9x  3x 

all  at  X  =  X 
o 


(22d) 


Note  that  the  bomdary  condition  given  by  equation  19  expresses  that  the 
shear  force  f  at  x  =  0  must  be  eqioal  to  F,  and  can  also  be  used  to 
calculate  F,  once  the  deformation  function  y(x,t)  has  been  determined. 
Note  also  that  F  can  be  obtained  by  ccxribining  Equations  7  and  8,  with 
the  condition  that  f  is  zero  at  the  downstream  end  of  the  ice  sheet,  and 
therefore  is  given  in  dimensionless  form  by 

^o  .2-  CO  or  t  2- 

F  =  I  (H  +  dx  +  J  (y  +  ~Il2^  (23) 

0  9t  Xq  8f 

or 


4F  = 


9X  3x 


X  =  0 


X  =  £ 


(23a) 


C.  Problan  Classification 

The  solution  of  the  governing  equation  will  depend  upon  the 
boundary  conditions  imposed,  vdiether  the  motion  can  be  assumed  to  be 
quasi-steady  or  not.  The  variolas  problans  are  schematically  represented 
in  Figure  3. 

The  problems  treated  in  the  following  sections  are 

1.  Quasi  steady  deformation  of  a  sooii-infinite  ice  sheet  for 

the  case  y  <  H  and  for 
o 

a.  F,  P  and  M  all  different  from  zero 

o 

b.  Only  F  different  fimi  zeie. 

2.  Semi-infinite  ice  sheet  in  quasi-steady  deformation, 

partially  emerged  or  submerged  (y^^  H),  with  F  ^  0  but  P  =  =  0. 
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3.  Quasi  steady  deformation  of  a  finite  ice  sheet  for  the 

case  y  <  H  and  with  F  #  0,  P  =  M  =0. 
o  '  o 

4.  Unsteady  deformation  of  a  semi-infinite  ice  sheet,  such 

that  y^  <  H  and  with  F  f  0,  P  =  =  0. 
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IV.  SOLUTIONS  OF  THE  QUASI-STEADY  EQUATION 


A.  Semi- Infinite  Ice  Sheet  in  Contact  with  Water 
Along  its  Vfliole  Length 

1.  General  Case  (F,  P  and  M  -/  0). 

o 

The  governing  equation  15  reduces  to 


.4-  .2- 

^+4(l)F^  +  4y  =  0  (24) 

dx  dx 


and  the  boundary  conditions  are 

as  X  ^  “ 

0 

X  =  0 

The  solution  of  equation  24  vdiich  satisfies  the  boundary  conditions  are: 
_X^x 

y  =  e  (k^  cos  A.2  ^  ^  ^2 


,n- 

y  and  ^  ^  0 
dx 


y  =  at  X  = 


,2- 


^  =  4  M  at 
-A  o 

dK 


with 

=  (1  - 


Frcan  either  equation  19  or  equation  23,  the  vertical  force  F  must  satisfy 
the  relation 


^2  =  (1  +  4>F)^ 
ko  = 


2M  +  (i)Fy 

O  "^O 


^1^2 


_  (l-2^F)y^  -2Mq 

2/1^ 


(26) 


The  bending  nonent  distribution  is  given  by 
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5  1  X  ,r 


(^l  )k^ 


COS  X 


■^  •  -  I 

sinA„x  \ 


_X,x  (_  _  M  (j)F  +  Yq/S 

=  e  |m  cos  X^x  + - :— - sin  X_  x 

I  o  2  X^Xg  2 


(27) 


The  location  x^,  of  the  raaxitnum  bending  mcanent  is  obtained  fron  ^ 
^  dx 

as 


x=x. 


=  0 


l+(j)F  (2(1)F-1)2M  +  y 


(2<))F+1)2M^  + 


(28) 


The  tensile  stress  at  either  the  top  or  bottom  of  the  ice 


sheet  at  the  section  x^  is  given  in  dimensionless  form  by 
0^  =  2  M^h  - 


(29) 


and  the  failijre  criterion  is 


=  Ob  (30) 

where 

a  =  a/E 

The  follomng  procedure  is  suggested  for  determining  the  forces 
F  and  P  exerted  by  an  ice  sheet  on  a  two  dimensional  sloping  structure. 

1.  The  geanetric  and  material  characteristics  of  the  ice 
sheet  (b,h,E,a^,  p^j^),  structure  angle  9,  and  the  friction  coefficient  u 
are  assumed  to  be  known.  Then  the  characteristic  length  L,  time  T,  and 
coefficient  (f),  are  known. 
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2.  Assume  a  value  for  and  therefore  =  yo/L. 

3.  Determine  the  corresponding  value  of  F  by  solving  equation 

26. 

4.  Determine  the  location  of  the  maximum  tensile  stress 
fran  equation  28. 

5.  Calculate  from  equation  27  and  frcra  equation  29. 

6.  Cotipare  with 

If  conputations  are  stopped 

If  repeat  steps  3  to  6  with  a  higher  value 

of  y  . 

•^o 


If  >  0|^  repeat  steps  3  to  6  with  a  analler  value  of  y^. 

It  should  be  noted  that  the  determination  of  F  frcm  equation  26,  requires 


that  be  estimated. 


P  by 


The  moment  M  is  related  to  the  hori2X3ntal  force 
o 


M  =  -  Pe 
o 

where  e  is  the  distance  fran  the  point  of  application  of  the  reaction 
force  R  to  the  neutral  axis  of  the  ice  sheet,  located  at  mid-thickness. 
When  the  ice  sheet  cones  in  contact  with  the  structure,  its  leading  edge 
will  be  sane\diat  crushed.  However,  the  eccentricity  e  can  be  approximated 
by  h/2,  then 


=  -  4.F  h/2 
M  =  -  (J>F  E/2 


(26a) 
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2.  Particular  Conditions  (F  /  0,  P  =  =  0). 

For  practical  application,  or  to  determine  a  starting  value  of 
y^,  it  may  be  sufficient  to  consider  that  the  ice  sheet  is  subjected 
only  to  a  vertical  force  F.  The  equations  derived  in  the  previous 
section  greatly  sinplify  to  read 


-  -X  - 

y  =  e  cos  X 


F  = 


o 


Xj  =  tan  ^1  =  ^ 
r.  - 


=  0.161  y 


o 


=  0.322  hy^ 


(31) 

(32) 

(33) 

(34) 

(35) 


and  therefore  the  force  at  failvire  (i.e, ,  vdien 
directly  in  terms  of  and  ii 

F  =  1.55  i 
^  h 

while  the  value  of  y^  at  failirre  is  given  by 


=  Oj^)  Ff  can  be  obtained 


(36) 


3. 


h 


(37) 


^^^re  ccranonly  encountered  dimensionless  forms  of  F„,  x  and  y  are 

I  I  O 

then  given  by 


Y  bh 
'  w 


^  w 


X 


^  =  0.60  (~)^  (-^)^ 


(38a) 


(38b) 
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It  is  WDrth  noting  that  according  to  equations  38,  the  force 
is  projxDrtional  to  the  bending  strength  and  inversely  proportional 
to  the  fourth  root  of  E,  only.  Therefore,  in  irodel  studies  of  ice 
forces  on  structures,  proper  scaling  of  the  bending  strength  is  rrore 
critical  than  that  of  the  elastic  modulus,  since  a  100  percent  error  on 
E  leads  to  an  error  on  F^  of  only  sixteen  percent. 

It  is  recalled  that  equations  38a  to  38c  are  valid  as  long  as 
the  free  end  of  the  ice  sheet  does  not  either  emerge  from  or  become  sub¬ 
merged  into  the  water,  i.e. ,  as  long  as 


(39) 


wiiere  a  is  equal  to  p./p  \^tien  the  sheet  is  lifted,  and  a  equals  (1  - 
p .  1  w 

— )  viien  the  sheet  is  pushed  down.  Therefore,  the  condition  of  appli¬ 


cation  of  equations  38a  to  38c,  can  be  expressed  from  equation  38c  as 


h  > 


3.20 


Y  E 


=  h 


max 


(40a) 


or 


E  3.20  . 

-  2  \ 

D  a  '  W 


(40b) 


3.  Exanple  of  Application 

Consider  an  ice  sheet  of  thickness  h  =  40  cm,  with  b  =  1  m, 
and  properties 

=  1  X  10®  Pa 

E  =  5  X  10®  Pa 

p^  =  920  kg/m®  p^^  =  1000  kg/m® 
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and  clirribing  up  an  inclined  structure  of  angle  0  =  30°.  We  will  assume 
a  friction  factor  between  the  ice  and  the  structure  of  y  =  0.15. 

It  can  be  verified  that  eq\xation  40  is  satisfied.  The  corresp¬ 
onding  characteristic  length  L  and  coefficient  (p  are,  respectively,  L  = 
10.21  m,  (j)  =  0.796. 

1st  Method:  If  the  effect  of  the  horizontal  force  P  and  mcraent  M  are 
-  o 

neglected,  then  equations  38a  to  38c  are  applicable  and  yield 

y  =  0.162  m 
°f 

=  8100  N 

Failinre  is  predicted  to  occur  at  a  distance  =  8.0  m  and  the  maximum 
horizontal  force  on  the  structure  is 

Pf  =  =  6450  N 

2nd  Method:  When  P  and  M  are  not  neglected,  but  under  the  assunption 
that  the  excentricity  e  is  equal  to  h/2,  i.e. ,  that  =  -  (J)F  h/2, 
the  problan  requires  the  simultaneous  solution  of  equations  26,  27,  and 
28,  varying  the  value  of  y^  until  the  bending  stress  given  by  equation 
29  is  fovind  equal  to  the  ice  bending  strength. 

With  the  nianerical  data  chosen  in  the  present  exanple,  the 
dimensionless  values  are  found  to  be 

y  =  1.635  10"^ 
o 

F  =  8.225  10“^ 

Xj  =  8.017  10"^ 


The  corresponding  dimensional  values  are 
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y  =  0.167  m 
F  =  8400  N 

=  8.19  m 
=  (t)F^  =  6700  N 

This  exanple  indicates  that  the  effect  of  P  and  on  the 
failvire  force  is  relatively  small,  of  the  order  of  5  percent  in  the 
present  case,  and  that  they  may  be  safely  neglected  and  still  obtain  a 
reliable  estimate  of  the  forces  exerted  on  the  structure. 


B.  Partially  Bnerged  or  Sutanerged  Sami-Infinite 
Ice  Sheet  (P  =  M  =0,  F  0) 

The  case  when  the  free  end  of  the  ice  sheet  becones  either 
fiilly  emerged  fran  or  fully  sutmerged  into  the  water,  i.e.,  for  y^  >  H, 
will  be  investigated  only  mder  the  assumption  of  quasi-steady  condition 
and  of  zero  horizontal  force  and  moment  at  the  upstream  ec^e  x  =  0. 

Under  these  assumptions  the  governing  eqmtions  16  and  15  reduce  to 

0  £  X  £  x^;  y^  >  H  (41a) 

x>x^-,y^<E  (41b) 

wiiere  x^  is  that  length  of  the  ice  sheet  \^ich  is  entirely  either  below 
or  above  water.  The  boimdary  conditions  are  given  by  equations  17  to  20 
and  22a  to  22d. 

The  solutions  of  equations  41a  and  41b  have  the  forms: 


9^1 


+  4  H  =  0 


9^yr 


+  4  y^  =  0 


3x 
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= 


s  ^  +  r  +  C„x^  +  r  X  +  C  , 
b  3  1  o 


(42a) 


X  <  X 

—  o 


rx' 


(42b) 


x'  =  X  -  X  >0 
o  — 


Note  that  equation  42b  satisfies  the  boundary  condition  given 
by  equation  20.  Prxxn  the  ramining  boundaiy  conditions,  the  coefficients 
k^,k2.  and  are  found  to  be 


=  H 


Cl  =  -  f  (2  + 


C2  =  0 


C3  =  I 


and  Xq  is  given  by  the  relationship 


_o_(l5^+2-3^-2^^  +1)=0 

5  ^6  o  3  o  o  o 

M 


(43a) 


or 


^  i  (d+x^)^  +  2(l+x^)  +  3) 
H 


O' 


o 


(43b) 
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The  force  F  exerted  at  the  free  end  of  the  ice  sheet  is  obtained 
from  equation  23,  for  quasi  steady  conditions,  as 

F  =  J  y  dx  =  H  dx  +  /_  ^2  dx 

0  X 

o 

1+x 

H  (44) 


The  ejqiressions  for  the  bending  nonent  M  obtained  fran  the 
dimensionless  form  of  equation  9  are  given  by 

—  H  —  —  _ 

M  =  75:  X  (1+x  -  x)  for  X  <  X 
^  o  o 

-  H  -  -  -  -  - 

M  =  (x  cos(x-x  )  +  sin(x-x  ))e 
^  o  o  o 

for  X  >  X 
o 


(45a) 

(45b) 


Failure  of  the  ice  sheet  occiars  \\hen  the  first  maximum  of  the 
bending  moment  satisfies 


M. 


(46) 


Examination  of  eqiiations  45a  and  45b  shows  that  the  maximum 

of  M  occTirs  at  a  distance  x^  less  than  x^,  i.e.,  in  the  fully  emerged 

or  submerged  part  of  the  ice  sheet,  \^ben  <  1,  but  greater  than  x^, 

for  x^  >  1.  For  illiistration  purposes  the  distribution  of  the  bending 

moment  is  shown  in  Figure  4  for  the  tvo  particular  values  x  =0.6 

and  X  =  1.4. 
o 

The  location  x^  of  the  maximum  bending  manent  corresponding  to 
(dM  /dx)-  -  =  0  is  then  given  by 

-1  -  -1  - 

x^  =  x^  +  tan  (-^)  =  x^  +  I  -  tan  (x^),  for  x^  <  1  (47a) 
1+Xo 

or 


f 


1+x 

o 


for  X  >  1 
o 


2 


(47b) 


0  0.5  1.0  1.5  2.1 


X 


Figure  4. 


Exanples  of_Bending  Moment  Distributio 
Slieet.  a)  x„  =  0.6;  b)  x  =  1.4 


o 


o 
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and  the  corresponding  ejqjressions  for  are 


/i  ^  -2 

a,  TT  /I  +  X 

M  =  —  =  —4/ - ^ 

f  2h  2 


ejq)  (-  ^  +  tan~^  x  )  for  x  <  1  (48a) 

^  4  o  o 


and 


a-1^  IT  1“*"^  2 

M  =  -^  =  I  (-^) 

^  2h  ^ 


for  X  >  1 
o 


(48b) 


Equations  48a  and  48b  are  plotted  as  2M^/H  versus  x^  on 


Figure  5.  When  x^  =  1  then 


Xj  =  =  1 


and 


M  -  H 
Mf  2 

The  results  of  the  analysis  for  the  case  x^>  1  deserve  further 

attention.  Coiribination  of  equations  44  and  48b  yields 


F  =  x^  H 


(49a) 


which  can  also  be  ejq^ressed  as 


- =  1  (49b) 

V  „bHx^ 

Furtherraore  by  eliminating  x^  between  equations  44  and  48b  one  obtains 
F  =  A  f  (50) 


which,  when  the  definition  for  L  is  introduced,  can  be  written 


F, 


Y  bh^ 

w 


H 

h 


(51a) 
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o 
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the  expression  for  can  be  reduced  to 


IT 

h  V  3 


%  h 


(51b) 


The  above  expressions  show  that,  vdien  x^/L  >  1>  both  F  and  are  independent 
of  the  elastic  modulus  E  of  the  ice,  and  therefore  of  the  ratio 
and  are  proportional  to  the  square  root  of  ^ly  the  displacanent 

y  will  be  function  of  E/q,  namely 

Oj  b 


!o  ,  8  hj  _b,  _b_  IT  ^  "TV 

h  g  V^E  \  h^  ^  ^  \h  ^E  >Y.h  2h 

3 


(51c) 


w 


The  condition  of  application  of  equations  51a  to  51c,  namely,  x^/L  >  1, 
can  also  be  written,  according  to  equation  48b,  as 


>  1 


(52a) 


hH 


or,  wiien  the  definition  for  L  is  introduced,  as 


a  1/2  0^  1/2 


,  “b  '  /  h 


(52b) 


that  is 


^  0.33  _  u 

-  ^2  y^E  min 


(53a) 


or 


E  0.33  °b  s 
0^  -  2r  V  h-' 

b  ct  '  w 


(53b) 


vdiere  a  =  H/h  and  can  be  ejq)ressed  in  terms  of  the  ice  specific  gravity. 

In  sunmary,  when  the  firee  end  fully  emerges  from  or  submerges 
into  the  water,  i.e. ,  for  h  >  h  as  defined  by  equation  40,  the  maximum 
forces  and  exerted  by  or  on  the  ice  can  be  determined  as  follows. 
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1.  Calculate  (j),  a,  L,  h  =  h/L,  H  =  H/L  =  ah,  and 

h  .  (equation  53a). 
mm  '  ^  ^ 

2.  If  h  >  h  .  ,  determine  x  fron  equation  49a,  x„  fron 

ram’  o  _  f 

equation  47a,  y  frcm  equation  43b,  and  F„  fron  equation  44,  Convert 

Oj^  I 

to  dimensional  quantities  x^,  and  F^, 

3.  If  h  <  h  .  ,  then  F„,  x-  and  y  are  obtained  directly  frcmi 

mm  f  f  o„ 

equations  51a  to  51c. 


4.  Calculate 


Exartple:  An  ice  sheet  of  thickness  h  =  15  can  is  being  piished  up  a  beach 
inclined  at  9  =  30°.  The  coefficient  of  friction  u  is  0.15,  and  the  ice 
properties  are 

“b  = 

E  =  2  X  10^  P 


Pi /p  “  0.92 
i'  w 


-  /  -  n  no  a  sm  9  +  y  CJOS  9 

1  w  ’  ^  cos  9  -  y  sm  9 


=  0.80 


^  4  Eh^  -  Q  QQ 

L  =Jo —  3,89  m 


H  =  ah  =  13.80  an 


_  3.20  ^ 

^max  2  Y  3  ^ 

a  'w 


''nmi  '  0-°2  ■" 


2,  Since  h  <  the  edge  of  the  ice  sheet  will  energe  frcm 
the  water  before  failure  occurs,  however,  h  >  h^^.^  and  the  equations  for 
x^  <  1  are  to  be  applied. 
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Equation  48a  yields  =  0.120 

Eqimtion  47a  yields  =  0.786 

Equation  43b  gives  y  /H  =  1.14 

Of 

_  _2 
Equation  44  gives  Ff  =  1.99  10 

The  dimensional  values  are 
Ff  =  2950  N/m 
y  =  15.7  cm 

Of 

Xf  =  3.06  m 


C.  Finite  Ice  Sheet  (P  =  =  0,  F  0) 

1.  Introduction.  We  have  examined  so  far  the  case  of  a  semi¬ 
infinite  ice  sheet  moving  against  a  sloping  structiire.  However,  in 
practical  cases  the  ice  sheet  is  seldom  semi-infinite.  Its  actual  size 
can  vary  frcan  that  of  a  river  ice  floe  of  relatively  small  dimensions 
striking  a  bridge  pier  or  other  structure,  to  that  of  a  lake  ice  cover 
of  very  large  dimension  being  pushed  up  a  beach  by  wind  shear.  We  will 
now  examine  the  case  of  a  finite  ice  sheet,  to  determine  the  length 
beyond  which  the  ice  sheet  can  be  considered  semi-infinite,  and  the 
equations  derived  in  the  previous  chapter  are  applicable. 

2.  Governing  Equation.  The  analysis  of  the  deformation  of  a 
finite  ice  sheet  of  finite  length  i  will  be  limited  to  the  quasi-steady 
case,  under  the  assuirption  that  the  lifting  edge  of  the  ice  sheet  never 
emerges  from  the  water  (i.e. ,  y^  £  h^)  and  is  subjected  to  only  a  lift 
force  F.  The  horizontal  force  P  =  0  and  therefore  =  0.  The  corresponding 
governing  equation  in  dimensionless  form  is 
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^  +  4y  =  0 
dx 


(54) 


The  edge  x  -  z  is  a  free  end  and  is  therefore  subjected  to  no 
moment  and  no  shear  force.  The  boundary  conditions  at  that  location  are 


i2-  .3- 

dy_dy  „  - 

dx  dx^ 


(55) 


At  the  lifting  edge  the  boundary  conditions  are 

,2- 

M  =  0  or  =0  at  X  =  0 


and 


(56) 


j3— 

d^ 

-3 

dx 


=  /  y  dx  =  4F 


x=0 


o 


(57) 


This  last  equation  results  frcm  the  combination  of  equations  23  and  23a 
with  equation  55. 

The  general  solution  of  equation  54  is  of  the  form 


y  =  e  cos  X  +  Kg  sin  X)  +  e^CKg  cos  x  +  sin  x)  (58) 


The  coefficients  K^,  Kg,  Kg  and  K^  are  obtained  in  terms  of  F  fran  the 
above  boundary  conditions  as 

K  =  -  sin(2t)  -  1 _  - 

cosh(2£)  +  cos(2£)  -  2 


^  4  cosh(2t)  +  cos(2£)  -  2  ^ 

jr  =  4  —  sin  (2£)  —  e _  ^ 

a  cosh(2£)  +  cos(2£)  -  2 

and  the  lifting  edge  displacanent ,  y^,  is  obtained  fixxn 
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y' 

•^o 


sinh(2&)  -  sin(2il) 
cosh(2Jl)  +  cos(2£)  -  2 


2F 


(60) 


A  sani-infinite  ice  sheet  subjected  to  the  same  lifting  force 
F  as  the  finite  ice  sheet  now  considered  has  a  di^lacanent  given  by 


'o 


=  2F 


therefore 


^o  _  sinh(2£)  -  sin(2£) _ 

cosh(2Jl)  +  cos(22,)  -  2  (61) 


Also  for  a  sani-infinite  ice  sheet  it  lias  been  shown  that  =  2F 
vtole 

increases,  tends  towards  2F  and  tend  to  zero. 

Table  1  presents  the  results  of  conputations  of  y^/y^,  K^/F 
(i  =  1,2, 3,4)  for  increasing  value  of  £.  It  can  be  seen  that  for  ^  = 
3.0,  and  are  practically  equal  to  y  and  2F  respectively,  and 
Kg,  K^  are  negligibly  small.  Therefore,  an  ice  flow  of  length 


Kg  and  K^  are  zero.  Equations  59  do  show  that  vlien  *- 


£  =  3.0L 


3.0  (|^) 
w 


3  1/4 


(62) 


can  be  considered  to  behave  as  a  soni- infinite  ice  sheet. 

^'^en  the  length  of  the  ice  floe  is  less  than  3.0L,  the  force 
exerted  on  the  structure  can  be  calculated  as  follows: 

1.  The  quantities  a^  =  K^/F  (i  =  1,2, 3,4)  are  calcxilated 
according  to  equations  (59). 

2.  The  location,  x^,  of  the  maximum  bending  nonent  M^,  is 
determined  fron  dM/dx  =  0  by  solving  numerically  the  equation 


(a^  -  a^)  siM^  -  (a^  +  a2)cosx^ 

3.  The  corresponding  value  of  M^/F  is  then  given  by 


(80) 


_f  -  1  d  y 
F  4F  dx^ 


X  =  x^ 
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T  -P 

=  cos  sinh  ^  (a^e^f  -  a^e  sin  (61) 

4.  The  maximum  bending  moment  at  failure  is  related  to  the 
bending  strength  of  the  ice  by 

\  =  2  h 

\diich  yields  the  ejqjression  for  the  failure  force 

Ff  =  _  (62) 

2h(Mj/F) 

5.  The  corresponding  displacement  of  the  lifting  edge  is  then 
obtained  as 

%  '  <®3> 


D.  Discussion 

Since  the  influence  of  the  horizontal  force  P  and  nonent  M  on 

o 

the  vertical  failure  force  appears  to  be  minor,  the  following  discussion 

is  limited  to  the  case  P  =  =  0,  and  is  expected  to  be  valid  when  P  and 

M  are  non  zero, 
o 

In  the  experimental  studies  and  modelling  of  forces  exerted 
by  monolithic  ice  sheets  on  inclined  structures,  in  particular  in  ice 
breakers  model  tests,  it  has  been  recognized  that  the  twD  parameters  E/a^ 
were  to  be  the  same  at  model  and  prototype  scales.  Frcan  the 
previous  analysis,  the  force  F  length  x  and  displacement  y  ,  for  the 

1  I  ^“F  ^ 

case  of  twe-dimensional  structure  can,  in  their  dinensionless^forms  F„/y  bh  , 

’  f''w  ’ 

x„/h  and  y  /h,  indeed  be  ej^ressed  as  functions  of  these  two  parameters 
only,  namely 


1.  For  E/a 


Y  h 

'  w 
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F  -1/4  a  3/4 

=  0-®  (f) 

Y  on  ^  w 

I  TIT 


(64a) 


if  =  (i-) 

b 


1/4  a  1/4 

'  ^17 


(65a) 


V  E  -1/2  °b  1/2 

-Tr-=1.79(~)  (-2^) 

h  \  \y 


(66a) 


with  also  the  relationship 


y 

2  °f 
IT  h 


O  T.  0.33  %  r^,  3.20  % 

2-  iWo 

a  ’  w  b  a  ^  w 


T 

_  0^  /L  ,  o. 

2  ^h  h  ^ 


(64b) 


T  r  1  X  h/ 

X-p  X  k  L  -  tan-l('-^  r) 
=  h  4  (h  1^ 


(65b) 


f  _  a 


I  [a+^)'‘  +  2(i+i)  +  3] 


(66b) 


h 

where  *  l  solution  of 


jl  +  (x  /L)' 

%  _  ,h.2  /  o 

E--  °‘(l)  V - i - 


+  tan 


T  X 


L._  1  ,E  .  "b  . 
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and  a  =  H/h 

No  ejqplicit  expression  similar  to  equation  (67 ) 
obtained  in  the  present  range  of  E/a^. 


For 


E/Ob 


0.33  “b 


Y  h 
'  w 


could  be 


3 


%  1/2 

V. 


(64c) 


If  =  /  1_ 

h  3a 


% 
'  w 


(65c) 


o  _ 


h 


8_ 

9a 


(^)(^)  +  0.358 

•tit 


'w 


0.50  a 


%  1/^  "b 


1/4 

'  TIT 


w 


(66c) 


vshile 


X 


o  _ 
h 


(70) 


and  also 


(71) 


For  future  use,  the 

were  plotted  versus  y*  =  /H  on  figure  6. 

An  infportant  resuft  is  that  ■when  the  ice  stiffness  S=  E/a, 

0  33  ^ 

IS  below  the  critical  value  S  =  - p-,  the  failure  force  F  and  the 

break  length  be(3ome  independent  of  S,  only  the  displacement  y^  romins 

dependent  upon  the  stiffness  F/oj^.  The  resilLt  has  inportant  implications. 

Consider  in  particular  the  dynamics  of  ice  breakers  and  the  problem  of 


quantities  F*  = 


X 

and 
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Figure  6.  Plot  of  F*  and  x  versus  y  /H 
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ice  breakers  model  tests.  There  exists  a  long  running  controversy  on  the 

proper  reduction  to  model  scale  of  the  mechanical  properties  of  ice.  The 

correct  scaling  of  the  ice  bending  strength  a  is  usually  achieved  by 

b 

using  doped  ice  (either  saline  or  urea  ice),  however  the  elastic  modulus  at 

scales  less  tlian  1/20  becomes  too  small  and  the  ratio  E/Oj^  in  the  model  is 

much  lower  than  that  in  the  prototype.  But  if  the  prototype  value  of  S  = 

E/ a,  is  less  than  its  critical  value  S  ,  then  the  model-ice  stiffness  could 
b  c 

be  lower  than  the  prototype  value  without  affecting  the  force  F^,  as  long  as 

the  model  ice  behaves  elastically  since  all  the  above  results  were  obtained 

assuming  that  elastic  theory  applied.  For  exaxiple,  if  for  sea-ice  the 

values  CT,  =  lO^ltPa,  p./p  =  0.91,  and  y  =  9810  N/m^  are 
b  1  w  w 

adopted^ since  ice  breakers  exert  a  downward  force  on  the  ice  the  corresp¬ 
onding  value  of  a  is  a  =  1  -  p . /p  =  0.09.  For  an  ice  thickness  h  of  2  m, 

F^  is  independent  of  E/a^  if  E/a^  is  less  than  2000;  \^en  E/Oj^  =  5000,  the 
force  F  will  become  independent  of  E/a,  for  ice  thicknesses  less  than  0.80  m. 

The  dimensionless  force  =  F„/y  bh  has  been  plotted  versus 
o  f'  w  ^ 

o  =  Oj^/y^h,  with  E/ a^  as  a  parameter,  on  Figure  7,  and  versus  E/oj^,  with 

as  a  parameter,  on  Figure  8.  In  these  figures  the  value  a  =  0.09  was 

selected. 


IV.  SOLUTION  OF  THE  UNSTEADY  EQUATION  (M^  =  P  =  1,  F  ^  0) 

Since  according  to  the  analysis  of  the  steady  form  of  equation 
15,  the  effect  of  and  P  on  the  horizontal  force  F  exerted  by  an  ice 
sheet  on  an  inclined  plane  is  relatively  minor,  only  the  case  =  P  =  0 
will  be  treated  here.  Furthermore  only  the  case  wben  failure  occurs 
prior  to  the  ccxtplete  eriergence  or  submergence  of  the  leading  edge  will 
be  considered.  The  gpveming  equation  is  then 


with  the  initial  conditions 


y(x,0)  =  0  all  X 


't 


=  0  all  X 


t  =  0 


and  the  boundary  conditions 


(73a) 

(73b) 
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_  an 

y  =  =  0  as  X  “ 

.-n 
ax 


ax 


=  0 


x=0 


^1 

at 


x=0 


„  -kt, 

"o '  ^ 

dt 


(74) 


(75a) 


(75b) 


where  k  is  sane  arbitrary  large  cx)ef f icient . 

The  form  of  the  boundary  condition  given  by  equation  75b  was 
chosen  so  that  it  can  be  coipatible  with  equation  73b,  and  that  the  ver¬ 
tical  velocity  of  the  leading  edge  reaches  a  practically  constant  value 
after  a  short  time. 

The  Laplace  transform  with  respect  to  time  of  equation  72  is 


with 


9x, 


^  +  4(1  +  s^)  y^  =  0 


(76) 


y^(x,s)  =L  ry(x,t)l  =  /” 

o 

The  boundary  conditions  becone 


y  = 

s 

Sn- 

— ^  =  0  as  X  “ 

3x“ 

(77) 

(78a) 

-2 

3X 

x=0  =  0 

y  (o,s)  =  y^  -  u  2 
®  so  O 

(78b) 

-  ■‘/  2 

Introducing  the  change  of  variable  x  =  xVl+s  into  equation 

S  ’ 

75  one  obtains 


+  4y 
•^s 


=  0 


(79) 
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vsiiich,  for  the  given  boundary  conditions,  has  the  solution 

u  k  -X 
„  _  o  _  s 

y„ - o -  e  cos  X 

®  s^(s+k)  s 


(80) 


According  to  eqioation  23  with  =  0,  the  force  F  is  given  by 

GO  2—  3— 

F(t)  =  /  (y+^)d5  =  l3^ 


o 


8t 


^  9x^  x=0 


the  Laplace  transform  F  of  F(t)  is  then 

s 

n  k(l+s^)3/4  u  k  (l+s^) 

•p  _  o  _ _  __  o 

s  r,.2,  „  2,  .  ,  ^  z 


2s  (s+k) 


2s  (s+k)  /l+s‘‘ 


(81) 


or 


u 
_c 

s  2 


F  =  -O 


1  A  -2^  4.  i+k 


s+k 


) 


(82) 


The  inverse  Laplace  transforms  of  the  various  functions  appearing  in  equa¬ 
tion  82  are 


.-1 


/tT 


r(i/4) 


t 

(f)  J 


-1/4 


(t)  = 


/tT 


r(i/4) 


G(t) 


^ere  (t)  is  the  Bessel  function  of  the  first  kind  of  order  i.  There¬ 
fore  the  force  F(t),  obtained  as  the  inverse  Laplace  transform  of  F  ,  is 

s  ^ 

given  by 


F(t)  =  J  [/G(X)dX  -  +  1±!£  G(,)g-k(t-T) 

2r(l/4)  ^  jo  k  '  ' 

Bessel  functions  of  the  first  kind  J  (x)  can  be  expressed  as  infinite 

A 

series  in  x,  namely 


dT(83) 
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c»  n 

J  fx')  =  F  (~1)  |,X\2n+A 

^  n!r(n+A+l)  ''2'' 


(84) 


therefore  the  f  met  ion  G(t)  can  be  expressed  as 


G(t)  =  (^) 


-  o=  C-l)^  t  2n-i 

J  1  q  (|)  2 

-1/4  n=0  n!r(n-f|)  ^ 


(85) 


When  this  expression  for  G(t)  is  introduced  into  eqioation  83,  the  various 
integrations  may  be  perfoimed  to  yield 


2u  /tT 


F(t)  =  —2—  E  - ^  v. 

r(l/4)  n=0  n!r(n+|) 


(-l)"^  ^ 


i (4n+l)(4n+3)  k(4n+l) 


k  p^O  p!(4n+2p+l)| 


(86) 


The  mcanent  distribution  in  the  ice  sheet  is  given  by 

8X 


the  Laplace  transform  of  v^ich  is 


M^(x,t)=  4-^  = 


2  -  £  2 

I  d  u„k/L+s  -X 


9x 


-2-5 -  e  ~  sin  x 

2s  (s+k)  ® 


(87) 


the  maximum  of  M  is  then  given  by 
s 


u 


kA+i 


M  „  =  0.161 
sf  2 


=  0.161  u 


s  (s+k) 


7:7  7 


1  ^  1+k  1 


k  s+k( 


(88) 


The  maximum  bending  mcxnent  is  obtained  as  the  Laplace  transform  of 


M.,  =  0.161  u 


00  iO 


J  (A)  dA 


o 


V  J  (t)  e 
k  o 


-k(t-T) 


dx  (89) 


When  is  replaced  by  its  infinite  series  ejqjansion 


y  /  -I  2n 

^  (~1^  xXs 


=  n=0  ^7^2  (f ) 


(n!)^ 


(90) 
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equation  89  becomes 


M„  =  0.161  u  ^  -2 

f  o  n=0  .  . .2 


n 


k‘ 


(niy 

li _ 

(2n-p)! 


|_A.  _  1)  , 

)2n+l  *^4 


1+k^  (2n)!  r2n  (-l)^(kt)^^  ^ 

LpSo 


-kt 


(90) 


which  can  also  be  written  as 


(-if  1. 


=  0.161  u_Io^ 


b 


n.  ,  2  X .  oo  1  p-2n-p 

1+k  (2n)!  .  .p  k^ _ 

^2  ^2n  p^2n+l  ^  p!  i 

”kt 

TAhen  the  expansion  of  e  is  used,  namely 

-kt  “  „  /'iJ\P 

e  =  E  (-l)’^-^IP- 


p=0 


(91) 


Failure  of  the  ice  sheet  occurs  when  the  maximum  bending 
mcment  satisfies  the  relationship 


= 


% 

2h 


(92) 


where  is  the  bending  strength  of  the  ice. 

Thus,  the  time  of  failure  t^  can  be  determined  frcm  equating 
the  right-hand  sides  of  equations  91  and  92.  Once  t^  is  determined  the 
force  exerted  by  the  ice  on  the  structure  is  given  by  equation  86  for  t 


tf. 
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VI.  EXPERIMEOTAL  VERIFICATION  AND  APPLICATION 

A.  Preliminary  Considerations 

The  foregoing  analysis  enable  to  predict  the  maximum  force,  F^, 
exerted  by  an  ice  sheet  on  an  inclined  structirre,  the  vertical  di^lacanent 
y^  at  failin:"e  of  the  edge  of  the  ice  sheet  and  the  length  x^  of  the 
brSken  floe  once  the  mechanical  properties,  and  E,  of  the  ice  are 
known.  Conversely,  the  strength  index  and  elastic  index  of  the  ice  can 
be  determined  from  lift-iip  or  push-down  ejqjeriments  in  vdiich  at  least 
twD  of  the  three  qxjantities  F^,  y^  and  x^  are  measured.  In  fact,  if 
these  three  parameters  are  measured  simultaneously,  then  each  pair  F^  - 
y  >  y  -  x„,  and  x„-F_  can  be  used  to  calculate  E  and  o  .  The  three 

O^  O^  III  1^ 

sets  of  E  and  deteirained  should  be  nearly  Identical,  within  the 
limits  of  ejqjerimental  accuracy,  if  the  assunptions  underlying  the 
analysis  are  satisfied,  in  particular  if  the  ice  behaves  like  an  elastic 
material . 

The  experimental  program  undertaken  at  the  Iowa  Institute  of 
Hydraulic  Research  to  verify  and  apply  the  previous  analyses  was  divided 
into  three  phases. 

Phase  1 . 

In  a  series  of  initial  ejq^eriments  the  bending  strength  and 
elastic  modulus  of  fresiiwater  ice  were  determined  and  conpared  with 
values  available  in  the  literature.  Freshwater  ice  under  sufficiently 
high  rate  of  strain  has  been  shown  to  behave  elastically  and  its  mechanical 
properties  are  well  documented. 

Phase  2. 

In  a  second  ejperimental  phase,  an  instrumented  inclined  plane 
of  variable  angle  of  inclination  was  built.  Forces  exerted  by  a  freshwater 
ice  sheet  pushed  at  constant  speed  against  this  structtire  were  measured 
and  coipared  with  the  theoretical  predictions  vdiere  the  values  of 
and  E  determined  in  the  initial  experiments  were  used. 
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Phase  3. 

The  results  of  experimental  phases  1  and  2  were  judged  satis¬ 
factory.  The  third  phase  was  concerned  with  the  determination  of  the 
strength  index  and  modulus  index  of  carbamide  (urea) -doped  ice  ■which  has 
been  proposed  by  Timco  (1979)  as  a  substitute  for  saline  ice  as  model  ice 
in  model  studies  of  ice  forces  on  struct'ures.  This  third  experimental 
phase  is  still  underway  and  only  preliminary  results  are  presented 
here. 

All  the  experiments  conducted  to  date  have  been  performed  in 
the  19  ft  X  3  ft  X  2  ft  (5.80  m  x  0.91  m  x  0.61  m)  ice  force  tank  of  the 
Low  Tetiperature  Flow  Facility  (LTFF)  of  IIHR.  The  ice  sheets  grown  in 
the  tank  were  made  free  frcm  the  tanlc  walls  by  placing  a  heating  tape  in 
a  narrow  sheet -metal  trough  mounted  at  water  le'vel  along  the  perimeter 
of  the  tank.  The  eiq^erimental  apparatus  and  procedures  specific  to  each 
series  of  tests  are  described  in  the  following  relevant  sections. 

Since  phases  1  and  3  above  are  concerned  with  the  experimental 
determination  of  the  mechanical  properties,  and  E,  of  ice,  the 
various  methods  of  data  analysis  used  to  that  purpose  and  based  on  the 
analytical  approach  developed  in  the  previoiis  chapter  are  first  presented 
in  the  following  section. 

B.  Methods  of  Calculations  of  and  E 

As  mentioned  in  Section  A  above,  the  mechanical  properties  of 

an  ice  sheet,  namely  the  bending  strength  index  and  the  elastic 

modiiLus  index  E,  can  theoretically  be  calculated  from  experiments  in 

which  the  free  floating  ice  sheet  is  subjected  to  a  pure  lift  or  push  at 

one  end,  and  when  two  of  the  three  parameters,  "the  displacement  at 

failiare  of  the  lifted  or  pushed  end,  F^,  the  force  at  failure  exerted  on 

the  ice,  and  x^,  the  breaking  length  of  the  ice  sheet  are  measured. 

It  is  assumed  that  the  following  data  are  known 

l :  initial  length  of  ice  sheet 

b :  width  of  ice  sheet 

h:  thickness  of  ice  sheet 

p  ^ :  ice  density 

p  ;  water  density 
w 
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therefore  the  quantity  H  can  be  calculated 


H  =  Pj^/P^  h  (lifting  case) 

H  =  (1  -  p  ./p,Jh  (pushing  case) 

X  w 


In  the  following  only  quasi-steady  conditions  are  considered, 
and  the  length  £  is  assumed  to  be  sufficient  for  the  ice  sheet  to  be 
considered  as  soni  infinite. 

The  results  of  the  preceding  analyses  with  =  P  =  0  are 
collected  below  for  future  reference 


/H<1 


also 


Sy  h 


Y  h  1/4 

=  0.68  bh 


O4 


Y  h  1/2  cv 
=  1.79  (— )  ~ 


E 


w 


Xf  =  0.60  h(^) 


1/4 


'W 


■p  _  2  -  ^  w*^^ 

^f  ^w^  ^f  ^o^  =  — T- 


(13) 

(38a) 

(38b) 

(38c) 

(93) 


Yo 

Of 


(43b) 


(44) 


(47a) 


(48a) 
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X  /L  >1 
o' 


1+  X  /L 
o' 


(47b) 


E 


H  h  /  o'  s 

L  L  ^  2  ^ 


(48b) 


that 


For  the  latter  case  (y  /H  >  1,  x  /L  >  1),  it  has  been  shown 


o 


o 


Y  bHx  „ 
'w  f 


=  1 


(49b) 


and  catbination  of  equations  13,  68b  and  69b  results  in 

Hx  /  X  ^ 

=  3y^ =  3aYw — 


(94) 


1.  The  y-F  method: 


The  deflection  y  and  force  F  are  measured.  The  unknowns 

I 

are  x^,  L,  E  and  Oj^,  and  x^'^when  y^  /H>1. 


a)  y  /H  <  1 :  The  characteristic  length  L  is  obtained  fron 


equation  93 


2F 


L  = 


Y  b  y 


(95) 


then  E  is  obtained  fron  equation  13  as 


E  = 


(S6) 


and  a,  and  x_p  are  given  by 

^  ^  F^2 

Y  bh^y  2 

'w  •'o^  2h 


(97) 
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ttF, 


2y  by 
w  Oj 


(38c) 


b)  y  /H  >  1;  the  ratio  x  /L  is  calculated  fron  equation  43b. 
o/  o 

Then  L  can  be  obtained  fron  equation  44,  E  from  equation  13,  fixxn 

equation  47a  or  49b,  and  0^  from  equation  69b  or  84,  depending  Aether 

X  /L  is  analler  or  greater  than  1,  respectively, 
o 


2.  The  x-y  method; 

The  deflection  and  length  x^  of  the  broken  floes  are 


measured. 


The  tinknowns  are  F^,L,  E  and 


and  X  if  y  /H  >  1 


Oj 


a)  yoj/H  <  1:  The  length  L  and  elastic  modulus  E  are  calculated 
from  38c,  fron  eqmtion  38b,  and  F^  from  equation  38a. 


b)  y^  /H  >  1:  x^/L  is  calculated  frc«n  equation  43b,  L  can 

then  be  determined  fron  equation  47a  if  x^/L  <  1,  or  equation  47b  if 

x^/L  >1.  F^  is  calculated  fron  equation  44,  and  E  fron  equation  13. 

Finally  is  calculated  fron  eqmtion  48a  if  x^/L  <  1,  or  frcxn  equation 

84  if  X  /L  >  1. 
o 


3.  The  x-F  method: 

The  floe  length  x^  and  force  F^  are  measured.  The  imknoms  are 
y^  ,  L,  E,  and  and  x^  vihen  applicable. 

^  The  quantity  F*  =  F^/f  calculated,  three  cases  are 

possible 


a.  F*  <  2/-n:  then  y  /H  must  be  less  than  1. 

Of 

determined  from  equation  38c,  from  equation  38a  and  y^ 
38b. 


L  and  E  are 

fron  equation 
f 


b.  2/tt  <  F*  <  1:  y  /H  must  be  greater  than  1,  but  x  /L 

o_p  o 

less  than  1.  Coribining  equations  44  and  47a  one  obtains  the  following 
equation  for  x^/L 


1 

2F* 


)  +  tt/4  — 


1  1 
2F* 

J 


which  yields  a  unique  solution  for  0  <  x^/L  <  1. 


(98) 
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Then  L  can  be  determined  frc*ii  equation  44,  y  from  equation 

43b,  E  from  equation  13,  and  fixan  equation  48a. 

c.  F*  =  1:  Then  x  /L  must  be  greater  than  or  equal  to  1, 

o 

and  /H  equal  to  or  greater  than  23/6.  The  only  quantity  ^diich  can 
be  determined  is  Oj^  from  equation  94.  The  elastic  modulus  E  is  undeter¬ 
minate  since  both  F^  and  are  independent  of  E. 

C.  Initial  Experiments 

1.  Experimental  Apparatus  and  Procedures: 

The  initial  series  of  tests  (Phase  1,  above)  consisted  in 
lifting  one  end  of  a  freshwater  ice  sheet  at  a  constant  speed  U^.  Each 
ice  sheet  tested  was  grown  in  the  LTFF  force  tank  at  a  room  tenperature  of 
-lO^C,  approximately.  Since  it  was  not  the  purpose  of  these  tests  to 
model  a  particular  ice,  the  ice  sheet  was  not  initially  seeded  but  left  to 
grow  naturally. 

The  instrumented  lifting  device  consisted  of  a  horizontal  angle 
iixon  panning  the  width  of  the  tank  and  su^ended  to  a  moment-insensitive 
dynamometer  by  two  vertical  bars.  The  dynamometer  was  attached  to  a 
vertical,  square-threaded  rod  passing  through  a  threaded  bevel  gear 
forming  one  part  of  a  right  angle  speed  reducer.  The  bevel  gear  was 
entrained  by  a  smaller  gear  mounted  on  the  shaft  of  a  variable  speed, 

D. C.  motor.  The  motor  speed  was  measured  by  a  light-interruptor  sensor 
set  over  the  perforated  edge  of  a  disk  mounted  on  the  motor  shaft. 

A  specially  designed  and  built  comter  comted  the  number  of  holes  sensed 
by  the  light-interruptor  and  converted  it  into  a  voltage  proportional  to 
the  vertical  displacanent  of  the  lifting  apparatus.  The  10-volt  full  scale 
of  the  counter  corresponded  to  5  ran  of  vertical  displacement.  The  counter 
reset  itself  to  zero  rtien  reaching  its  full  scale,  so  that  the  vertical 
displacement  was  obtained  as  multiple  of  Snnn  incranents  and  portions. 

The  vertical  force  F  exerted  on  the  lifting  devices  was  measured  by  a  com¬ 
bination  of  a  load  cell  (Statham  T-'/UL-d)  and  a  transducer  (Statham  iWC-2) 
mounted  on  the  dynamometer.  Both  the  vertical  displacanent  y^  and  the 
force  F  were  continuously  recorded  on  a  Beckman  multiple  channels  recorder. 
A  sketch  of  the  experimental  apparatiis  is  shown  in  Figure  9. 

Failure  of  the  ice  sheet  was  indicated  by  a  sudden  drop  in 
the  force  record.  Maximum  force,  F^,  and  corresponding  displacement. 


TO  RECORDER 
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,  were  measured  fron  the  records.  The  length  of  the  broken  floe,  x^, 
was  also  measured.  From  these  data,  the  characteristics  and  E  of  the 
ice  were  evaliiated  by  the  three  iiBthods  discussed  in  the  preceding  section. 

2.  Experimental  Results 

Six  freshwater  ice  sheets  of  various  thicknesses  were  tested. 

Since  all  the  ice  sheets  were  grown  in  the  same  manner,  their  mechanical 
properties  ought  to  be  the  same  within  the  limits  of  e^^rimental  acciiracy 
and  validity,  of  the  analysis. 

The  experimental  data  and  results  of  the  analyses  are  presented 

in  Table  2.  Since  all  the  sheets  were  relatively  thin  and  stiff,  failin:*e 

always  occurred  after  the  lifted  end  had  fully  emerged  from  the  water. 

All  sheets  had  a  width  b  of  0.80  m  and  an  initial  length  £  of  4.40  m. 

From  the  results  listed  in  Table  2,  it  can  be  seen  that  the 

calctilated  values  for  a ,  given  by  the  three  methods  are  very  consistent 

b  5  2 

and  renain  within  ten  percent  of  an  average  value  of  10.9  x  10  N/m  . 

On  the  other  hand,  the  x-F  method  gives  usually  larger  values  of  E  than 

either  the  y-F  or  the  x-y  method,  both  of  ^^ch  predict  E  values  nearly 

equal.  The  y-F  method  predicts  an  average  value  for  E  of  (1.74+0.35) 

X  10^  N/m^,  the  x-y  method  a  value  of  (2.06  +  0.37)  x  10^  N/m^,  and  the 

9  2 

x-F  method  a  value  of  (4.0  +  1.9)  x  10  N/m  .  Since  E  is  proportional 
4 

to  L  ,  relatively  small  discrepancies  in  the  calculated  value  of  L  lead 
to  large  variations  in  E.  As  can  be  seen  from  figure  6,  in  the  neighborhood 
of  F*  =  1,  small  changes  in  F*  correspond  to  large  variations  in  x^/L 
and  Fq/H,  and  therefore  in  L;  for  exanple,  as  F*  varies  from  0.95  to 
1.0,  a  change  of  only  five  percent,  x^/L  varies  from  0.6  to  1.0,  a  forty 
percent  change,  and  y^/H  varies  fron  2.1  to  3,8,  an  eighty  percent 
change.  To  illustrate  f\n:i;her  the  above  remarks,  let  las  consider  experiment 
F-1.  With  the  given  values  of  h,  F  and  Xf ,  F*  is  calculated  to  be 

9 

0.989,  and  the  corresponding  values  of  L  and  E  are  0.77  m  and  4.6  x  10 
2 

N/m  ,  respectively.  However,  x^  can  be  measured  at  best  within  one 

centimeter,  since  the  break  in  the  ice  sheet  is  never  perfectly  straight, 

and  x^  represents  an  average  floe  length  over  the  width  of  the  ice 

sheet.  F*  could  have  been  taken  equal  to  1.0,  just  as  well;  the  corresponding 

9  2 

values  for  L  and  E  would  have  been,  L  =  0.69  m  and  E  =  3.0  x  10  N/m  . 


Table  2.  Calculation  of  Ice  Mechanical  Properties  -  Initial  Experiments 
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In  conclusion  to  the  above  disciissions ,  it  is  suggested  that  the  x-F 
method  be  used  to  determine  the  elastic  moduliis  index  E  only  ■vdien  the 
quantity  F*  is  less  than  0.85  to  0.90.  For  values  of  F*  greater  than 
0.85,  the  y-x  and  y-F  methods  are  to  be  preferred,  even  though  it  is 
likely  that  x^  and  F^  are  more  easily  measured  than  y^^.  All  three 
methods  can  be  used  to  obtain  reliable  estimates  of  the  strength  index 

%• 


D.  Experiments  with  Sloping  Planes 

In  the  second  experimental  phase  of  the  study,  a  sloping 
structure  was  mounted  in  the  tank.  Its  slope  coifLd  be  varied  from  10° 
to  80°.  This  structure,  constructed  of  aluminum  angles,  spanned  the 
entire  width  of  the  tank  and  was  instrumented  with  two  moment  insensitive 
dynamcraeters  to  measure  the  horizontal  and  vertical  force,  F  and  P 
exerted  by  an  ice  sheet  pushed  at  constant  speed  against  it.  A  si^etch 
and  photograph  of  the  instrumented  sloping  plane  are  shown  in  Figures  10 
and  11,  respectively. 

Five  ejqjeriments  were  conducted,  three  with  the  plane  Inclined 
at  0  =  37°,  and  two  at  9  =  45°.  In  each  test,  continuous  records  of  the 
forces  F  and  P  were  obtained  fron  which  the  failure  forces  F^  and  P^, 
and  the  time  of  failure  t^  were  measured.  In  addition  the  length  of  the 
broken  floe  was  also  measured. 

Since  the  ice  sheets  used  in  these  experiments  were  grown  in 
the  same  way  as  those  tested  in  the  first  experimental  phase,  the  value  of 
adopted  in  the  prediction  of  F^  and  P^  was  the  average  of  all  values 
listed  in  Table  1,  vhile  that  of  E  was  the  average  of  the  values  of  E  in 
Table  1  predicted  by  the  Y-F  and  X-Y  method,  namely 


=  10.9  X  10^  N/m^ 
E  =  1.9  X  10^  N/m^ 


The  theory  developed  earlier  for  the  conditions  P  =  0,  M  ==  0, 
F  ^  0,  was  applied  to  predict  the  value  of  F^,  the  corresponding  value 
of  P^  was  then  calculated  as 


P„  =  F„  tan  0 
f  f 
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Apparatus  in  operation 


Close-up  view 


Figure  11.  Photographs  of  Eb^rimental  Sloping  Plane 


56 


that  is  that  the  friction  factor  y  between  ice  and  structure  was  assumed 
to  be  zero.  This  approximation  was  based  on  the  fact  that,  by  construction 
of  the  experimental  structure,  the  ice  sheet  was  supported  by  only  five 
smooth,  narrow  aluminum  angles,  as  can  be  observed  on  Figure  11.  Also 
the  angle  0  coiild  be  measured  at  best  within  an  acciiracy  of  one  degree, 
and  therefore  tan  0,  for  0  =  37*^  or  45°,  was  known  within  an  accuracy 
of  +  4  percent.  Finally  repeated  calibration  of  the  force  measuring 
systans  indicated  an  acciiracy  in  the  measured  values  of  F^  and 
of  approximately  5  percent. 

The  ejqjerimental  data  and  the  predicted  values  of  F^,  P^,  and 
Xj  are  listed  in  Table  3.  All  the  ice  sheets  were  thin  enough  for  the 
riding  edge  to  emerge  fron  the  water  since,  the  minimum  thickness  for 
the  sheet  not  to  anerge  is 


% 


P  2 

h  =  3.2  (— )  =  0.24  m 

1  '  W 


vdiich  is  one  order  of  magnitude  larger  than  the  thickness  of  the  experimental 
ice  sheets. 

From  the  results  of  Table  3,  it  can  be  seen  that  the  predicted 

values  of  F^  are  in  very  good  agreanent  with  its  measured  values  with  a 

maximum  difference  of  ten  percent  for  the  last  and  thickest  ice  sheet. 

The  calculated  values  of  P  aare  in  good  agreanent  with  the  measured 
o  ^ 

ones  for  0  =  45  ,  but  are  twenty  to  twenty  five  percent  smaller  than  the 
measured  values  for  0  =  37°.  Besides  experimental  errors,  it  is  likely 
that  part  of  the  discrepancy  is  due  to  the  neglect  of  friction.  For  0  = 

37°,  then  <j)  =  0.75  if  y  is  assumed  to  be  zero.  The  average  measured 
value  of  (j>  for  0  =  37°  is  0.97  corresponding  to  a  value  of  y  =  0.13. 

With  the  latter  value  of  y,  the  calculated  P^  would  fall  within  five 
percent  of  the  measured  one. 


E.  Experimental  Determination  of  and  E  for  Urea-Doped  Ice 

Tlmco  (1980)  proposed  that  urea-doped  ice  oe  vised  instead  of 
saline  ice  as  model  ice  in  laboratory  modeling  of  ice-structure  interac¬ 
tion  such  as  ice  brealcers  model  studies.  The  results  of  Ills  ej^riments 
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on  urea-doped  ice  indicate  tliat  the  ratio  E/a  remains  equal  to  or  larger 

b 

than  2000  for  a  flexural  stren^h  as  low  as  20  kpa,  for  ice  grown  from 
a  1.3%  solution  of  carbamide  (urea)  in  water.  These  values  of  the  mechani¬ 
cal  properties  of  the  model  ice  would  allow  a  model  scale  up  to  1  =50, 
while  the  maximum  value  of  X  is  of  the  order  of  20  wben  saline  ice  is  used 
as  model  material.  In  Timco's  ejqjeriments,  the  strain  modulus  E  was 
measured  by  the  plate  method,  and  the  flexural  strength  was  measured 
by  the  cantilever  beams  method,  with  the  beams  loaded  downward. 

It  was  then  decided  to  determine  the  mechanical  properties  of 
urea-doped  ice  by  the  method  of  the  floating  ice  sheet  described  in 
section  VI -C  above  and  previously  applied  to  freshwater  ice  sheets. 

A  first  series  of  ejqjeriments  was  conducted  with  the  lift-piish  apparatus 
used  in  the  freshwater  ice  experiments.  This  apparatus  permitted  to 
achieve  a  maximum  vertical  speed  of  approximately  3.5  ran/ sec.  Because  of 
concern  that  at  large  urea  concentration  the  rate  of  deformation  might  not 
be  high  enough  to  insinre  elastic  deformation  of  the  ice,  the  apparatias 
was  later  modified  by  eliminating  the  right -angle  speed  reducer  and  connec¬ 
ting  the  motor  shaft  directly  to  the  threaded  rod  supporting  the 
instrumented  lift -push  device.  The  maximum  vertical  speed  was  then  in¬ 
creased  to  about  30  ran/ sec.  A  photograph  of  the  direct  drive  lift-piish 
apparatus  is  shown  in  Figure  12, 

The  ejqjerimental  program  was  designed  to  investigate  the  effect 
of  initial  urea  concentration  in  the  water  bath,  test  tenperatxire ,  and 
load  direction  (piish  or  lift).  The  tests  perfomed  so  far  have  all  been 
conducted  in  the  lift  mode,  and  only  three  urea  concentrations  have  been 
tested,  namely  0.3,  0.6  and  0.75%.  Therefore  the  results  presented 
in  Table  4  are  to  be  considered  as  preliminary.  The  general  experimental 
procedure  was  as  follows. 

1.  The  air  tQiperatxu^  was  brou^t  to  about  -7.5°C  and  the 
water  bath,  except  for  the  first  sheet  tested,  was  constantly  circulated 
in  order  to  achieve  a  honogeneous  taiperature  tliroughout.  When  the 
initial  ice  layer  formed  naturally  at  the  water  surface,  the  layer  was 
slcirrmed  off  and  the  bath  was  wet  seeded  by  spraying  a  fine  mist  above 
the  tank.  The  resulting  ice  crystals  formed  in  the  air  deposited  at  the 
surface  of  the  water  and  initiated  the  test  ice  sheet.  The  ice  sheet 
Was  grown  for  twenty  foiir  to  thirty  six  hoinr’s  to  reach  a  thickness  of 
several  centimeters. 
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Front  View 


Back  View 


Figure  12.  Photograph  of  Direct  Drive  Lift-Push  Apparatus 
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2.  Once  the  ice  sheet  had  reached  the  desired  thickness  it 
was  either  tested  irnnediately  at  a  roan  toiperatijre  of  -7. 5*^0,  or  the 
refrigerated  systan  was  turned  off  and  the  roan  temperatiare  let  to  rise 
slowly  to  about  -1°C  at  vdiich  time  the  tests  were  conducted.  The  warm 
up  tine  varied  fron  twelve  to  eighteen  hours. 

No  thin  sections  of  the  ice  were  made  to  examine  its  structural 
features;  however,  it  was  observed  that  the  urea  ice  was  quite  unhono- 
geneous  with  a  hard,  thin  top  layer  veiy  similar  to  freshwater  ice, 
and  a  thicker,  mushy  botton  layer  similar  to  saline  ice.  Overall,  urea 
ice  appeared  to  be  more  "poroiis”  than  saline  ice;  ^hen  remeved  fron  the 
melt,  an  ice  sanple  weuld  be  drained  very  rapidly  of  its  liquid  constituent. 
Also,  freeboard  observations  indicated  that  urea  ice  is  denser  than  both 
freshwater  and  saline  ice  in  agreement  with  Timco’s  observations.  This 
latter  feature  became  apparent  for  \nrea  concentration  of  0.5%  and  larger. 

The  ejqjerimental  records  of  lift-force  versus  time  often 
showed  the  force  increasing  then  leveling  off  with  no  sudden  drop,  wiiich, 
in  the  case  of  freshwater  ice,  was  indicative  of  failure.  However, 
a  line  of  fracture  across  the  ice  sheet  could  be  observed  to  form  approxi¬ 
mately  at  the  time  the  lift  force  started  to  level  off.  Because  of  this 
feature,  the  displacement  y^  of  the  free  end  at  failure  could  not  be 
determined  with  certainty,  ft  was  then  decided  that  and  E  wDuld  be 
determined  by  the  X-F  method,  taking  as  the  maximum  force  recorded 
and  Xj  the  length  of  the  broken  floe. 

The  experimental  data  and  the  results  of  the  coiputations  are 
shown  in  Table  4.  It  can  be  observed  that  oj^  decreases  slowly  with 
increasing  tnrea  concentration,  is  only  slightly  affected  by  the  air 
tonpemture,  and  that  its  values  for  noninally  identical  experimental  con¬ 
ditions  remain  reasonably  constant,  with  maximum  deviation  of  about  50%. 

On  the  other  hand,  the  calcxiLated  values  of  E  can  vary  widely  between 
experiments  conducted  under  the  same  noninal  conditions.  It  can  also  be 
noted  that  for  the  same  urea  concentration  and  air  taTperat\ire ,  the 
thinner  ice  sheets  yield  larger  values  of  both  Oj^  and  E  (for  exanple, 
ice  sheets  and  6  or  7).  It  can  be  conjectured  that  most  of  the  strength 

and  stiffness  of  the  ice  sheet  arises  fron  the  top  layer  vdiich  occupies 
a  larger  percentage  of  the  total  ice  thickness  for  thin  ice  sheets. 

As  mentioned  previously,  the  experiments  with  urea  doped  ice  are 
still  underway,  and  mere  coiprehensive  results  will  be  available  \dien  the 
experimental  program  is  coipleted. 


Table  4.  Results  of  Experiments  with  Urea-Doped  Ice 
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VII.  SUMMARY  AND  CONCLUSIONS 

Under  the  assunption  that  ice  is  a  hotnogeneous ,  elastic  material, 

analytical  expression  for  the  hori2iontal  and  vertical  forces  exerted  by 

an  ice  sheet  on  a  two-dimensional  sloping  plane  have  been  derived.  It 

has  been  shown  that  an  ice  sheet  of  length  three  times  or  mere  a  charac- 

3  ^ 

teristic  length  L,  defined  as  (Eh  /Sy  )^,  can  be  considered  as  seni-infinite. 

Also  the  analytical  results  indicate  that  the  horizontal  force  and  bending 

moment  at  the  frontal  edge  of  the  ice  sheet  have  little  influence  on  the 

deformation  of  the  ice  sheet.  The  ice  sheet  can  then  be  treated,  for 

practical  application,  as  subjected  only  to  a  vertical  force.  The  resulting 

analysis  and  final  expressions  for  the  maximum  force  exerted  on  the  sloping 

plane  are  considerably  sitiplified  without  significant  loss  in  accuracy. 

Two  critical  ice  thiclmesses  have  been  identified,  h  .  and  h 

’  mxn  max 

They  have  been  expressed  in  terms  of  the  mechanical  properties  and  E  of 
the  ice,  and  its  specific  gravity,  and  have  been  shown  to  depend  on  the 
direction  of  application,  upward  or  downward,  of  the  vertical  force  F  exerted 
at  the  frontal  edge  of  the  ice  sheet.  IVhen  the  ice  thickness  is  larger  than 
Hiax’  failure  force  F„  is  proportional  to  a^,  inversely  proportional  to 

E4,  and  proportional  to  ^Vhen  the  ice  thickness  is  less  than  l\nin>  ^-e 

h  312,  ^ 

is  proportional  to  and  to  h  '  ,  but  independent  of  E.  This  latter  result 

has  the  irrportant  consequence  that  in  model  studies  of  ice  forces  on 

inclined  structures,  ^dlere  the  condition  h  <  h  .  is  satisfied,  only  the 

bending  strength  of  the  ice  need  be  properly  scaled,  v^hile  the  elastic 

modulus  of  the  model  ice  may  take  any  value,  as  long  as  the  model  ice 

behaves  as  an  elastic  material. 

The  analytical  expressions  have  been  verified  satisfactorily 

in  experiments  conducted  with  freshwater  ice  sheets.  An  experimental 

pixogram  has  been  initiated,  ^^here  the  analytical  expressions  are  used  in 

the  determination  of  the  mechanical  properties  of  urea-doped  ice.  Only 

preliminary  results  of  this  ejgDerimental  program,  which  is  still  undeiway, 

have  been  presented. 

It  remins  to  be  verified,  either  analytically  or  experimentally, 
that  the  main  results  obtained  here  ramin  valid  for  the  case  of  ice  forces  on 
tliree -dimensional  structures  such  as  sloping  planes  of  finite  span  and  conical 
structures.  In  particular  if  there  exists  a  critical  ice  thickness  below 
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’wiiich  the  maximuim  ice  force  is  independent  of  the  ice  elastic  modulxis, 
model  studies  could  be  perfoimed  without  concern  about  the  elastic  modulus 
of  the  model  ice. 
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